The photon polarization in B 



X7 in the standard model 



o 
o 

(N 
o 

Q 



> : 
o\ . 

o ■ 

<n : 

t-H . 

-xf ■ 

o ; 
^ : 

Oh! 

1 

Oh! 



X 



Benjamin Grinstein, 1 Yuval Grossman, 2 ' 3 ' 4 Zoltan Ligeti, 5 and Dan Pirjol 6 

1 Department of Physics, University of California at San Diego, La Jolla, CA 92093 
2 Department of Physics, Technion-Israel Institute of Technology, Technion City, 32000 Haifa, Israel 
s Physics Department, Boston University, Boston, MA 02215 
* Jefferson Laboratory of Physics, Harvard University, Cambridge, MA 02138 
5 Ernest Orlando Lawrence Berkeley National Laboratory, University of California, Berkeley, CA 94720 
6 Center for Theoretical Physics, Massachusetts Institute for Technology, Cambridge, MA 02139 

The standard model prediction for the B° — > X Sj( j7 decay amplitude with a right-handed photon 
is believed to be tiny, suppressed by nig^/mb, compared to the amplitude with a left-handed photon. 
We show that this suppression is fictitious: in inclusive decays, the ratio of these two amplitudes is 
only suppressed by g 3 /(4ir), and in exclusive decays by Aqcd/ti6- The suppression is not stronger 
in B° — » Xd~i decays than it is in B° — > X s j. We estimate that the time dependent CP asymmetries 
in B — > K*^j, p7, ifs7r°7, and 7r + 7r _ 7 are of order 0.1 and that they have significant uncertainties. 



I. INTRODUCTION 

The standard model (SM) predicts that photons are 
mainly left-handed in b — > (77 (q = s, d) decay (and right- 
handed in b — > 57). We define the ratio 



r a e 



Ar = A(B 



fqlR) 



A(B -» f qlL ) 



(1) 



where is a weak phase and 5 is a strong phase. It 
is usually stated that r q = m q /mb < 1 in the SM [1]. 
(Throughout this paper B refers to B° or B~ that con- 
tain a b quark, and r, 4> and S depend on the final state, 
/.) New physics can modify this prediction, and there- 
fore several methods have been proposed to measure the 
photon helicity [1, 2]. 

In B — > /7, where / is a CP eigenstate, since 7 l and 
7_r cannot interfere, the time dependent CP asymmetry 



r[g°(t)-+/7]-r[B (f)^/ 7 ] 
r[B°(t)->/ 7 ] + r[B°(t)->/ 7 ] 

= S/ 7 sin(Am t) — Cj 1 cos( Am t) , 



(2) 



is sensitive to r. In the SM, (f> s and C/ s7 are suppressed 
by \{V ub V us ) / (V tb V ts )\, and to first order in r s (< 1) 



Sf sl — —2r s cos 5 S sin 2(3 . 



(3) 



The first measurements of such CP asymmetries were 
carried out recently, 
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BELLE [4] , 
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yielding a world average Sk*j = —0.28 ± 0.45. At a 
super- B- factory the statistical error with 50ab _1 data 
is estimated to be 5(S K *-y) = 0.04 [5]. The Belle Col- 
laboration also measured the CP asymmetry S'^- s7r o 7 = 



-0.58 



+0.46 
-0.38 



± 0.11, integrating over the invariant mass 



range 0.6 GeV < m Ks7T o < 1.8 GeV [6]. It will also 
be possible to measure this CP asymmetry in B — > 
7r+7r - 7 [7], and maybe even with additional pions [8]. 

The purpose of this paper is to study the SM predic- 
tion for r. (For earlier attempts to go beyond the naive 
estimate, see Refs. [9, 10].) We find that r is only sup- 
pressed by <7 s /(47t) in inclusive b — > A7 decay, and by 
Aqcd/wi6 in exclusive B — > K*j and p7 decay. 

To understand the origin of such effects, recall that the 
effective Hamiltonian for b — > 57 is [11] 



H Q s — — - 



AG} 



-IVMYjCiitiOiin). (5) 
For our discussion the operators directly relevant are 



2 - (c^P L b)(s^P L c), 



Oj = 



I67P 



(6) 



where Pr.l = (1 ± 7s)/2, and we neglect the m s PL part 
of O7 hereafter. At the parton level, as long as b — > sj is 
a two-body decay (either from the leading contribution 
of O7 or subleading virtual contributions from the 
left-handed s quark is back-to-back to a photon. Then 
the two-body kinematics implies that only 7^ is allowed. 
This argument does not apply to multi-body final states, 
such as b — * 57 + gluons. 

The TObPfj part of the leading operator O7 contributes 
only to Al to all orders in the strong interaction. To 
prove this, note that the electromagnetic tensor for jl,r 
is Ffc R = \{F^ ± iF^), where F^ = \e llvpX FP x , and 
O7 can be written in terms of mbF^ v . Thus, independent 
of hadronic physics, the photon from O7 is left-handed. 
This argument only applies for O7. Indeed, we find by 
explicit calculation that other operators produce right- 
handed photons once QCD corrections are included. 
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FIG. 1: Leading contribution to B — > X s ^r. A second dia- 
gram with the gluon and the photon exchanged is not shown. 



II. INCLUSIVE B -f X s y AND B -f Jf d 7 

In this section we estimate r from an inclusive calcu- 
lation. The result can only be trusted if several hadronic 
final states are allowed to contribute, and r for specific fi- 
nal states cannot be obtained from this calculation model 
independently. 

The leading contribution to the inclusive B — > X s -yn 
rate is of order a s . It arises from bremsstrahlung con- 
tributions to the matrix elements of operators other that 
Oj. The numerically dominant contribution comes from 
O2 shown in Fig. 1. The corresponding amplitude was 
calculated in Ref. [12, 13]. We find that it yields equal 
rates for left- and right-handed photons at order a s , at 
any point in the b — > sgj Dalitz plot. 

Because of the complicated m c -dcpendence of the dou- 
ble differential rate, dr^ 6 ™^ /dE^dEg, we integrate over 
Eg and E 1 numerically. To reduce the large sensitivity to 
the scale of a s , we include the known order a 2 Po contri- 
bution to T^ em ^ [14]. (At this order the equality of the 
decay rates to left- and right-handed photons is violated 
by less than 1%, and can be neglected.) 

Using the "effective" Wilson coefficients at leading or- 
der [15], C2(mb) = 1.1 and CV(mb) = —0.31, a s (mb) — 
0.22 with to;, = 4.8 GeV, and m c — 1.4 GeV, we obtain 
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To 



0.025, 
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This result corresponds to integrating the numerator over 
x = 2E 1 /m b > 0.75 (that is, roughly, E 1 > 1.8 GeV). 
This result includes also the 0{a 2 8 [3$) correction, which 
is sizable, indicating that the relevant scale of a s may be 
well below mt; without including it the result in Eq. (7) 
would be 0.015. Thus, we find at lowest order in g s 



can chose the phase conventions independently. These 
amplitudes can still contribute to the same hadronic fi- 
nal states and interfere once hadronization effects are in- 
cluded. Thus, the relevant phase for any final state is de- 
termined by the hadronization processes and cannot be 
extracted from the inclusive calculation. Comparing the 
absorptive and dispersive parts of the inclusive result, wc 
find cos 2 S s — 0.3 with small variation over 0.75 < x < 1. 
The only conclusion we can draw here is that we expect 
the strong phase to be generically large. 

Next we discuss inclusive B — > X d j. In Eq. (5), op- 
erators multiplying the suppressed CKM factor V u bV* s 
were neglected. The analogous terms are important in 
b — > dj, since in this decay the u quark loop is not CKM 
suppressed compared to the c and t loops. The b — > dj 
effective Hamiltonian is of the form 



AC 

* i=l 



(9) 



AG F 



where prime denotes that in the standard operator basis 
in Eq. (6) the s quark is replaced by d, and the u subscript 
means that c and c are replaced by u and u; for example, 



'2,u 

It is useful to define 



0' 2 , u = (^PLb) (d^P L u) 



i{<t>d+$d) 



i8 3 , VubV* d u i5 « 
s T/ T/* d ■ 



(10) 



(11) 



where the first [second] term comes from the contribution 
of the top [bottom] line in Eq. (9) to Ar. The first term is 
the same as in b — > s decay (recall that we neglected m s ) . 

For the u quark loop, taking the limit of the quark 
mass to zero, the calculation of Ref. [14] simplifies con- 
siderably, and we can obtain analytic results 



1 dr. 
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— ln[x(l - x)] + — In x 
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x>0.75 
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where a s is evaluated at the scale nib in the MS scheme. 
(8) For the difference of the rates to -fa and 7^ we obtain 



This value of (r s ) decreases only slowly with a stiffcr cut 1 d[T ( ^ R) - T [ £ L) ] _ C\ a 2 s (3 x 2 (3 - 3x + 4x\nx 



on x. 

The value of cos<5 s is physical, as it enters S/ s7 in 
Eq. (3). Yet it cannot be estimated from the inclusive 
calculation. The reason is that the dominant contribu- 
tion to A R comes from the b sjg amplitude generated 
by O2, while to from the b — > 37 decay generated 
by O7. These are different final states, for which one 



dx 



Cf 27tt 2 



(13) 

This difference integrates to zero, and it gives a slight 7l 
enhancement near x=l. We obtain for x > 0.75 



^(u,brem) 
22 



To 



~ 0.030. 



(14) 
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The absolute values of the amplitudes corresponding 
to the c and u loops [the squares of which yield Eqs. (7) 
and (14)] are comparable to each other. In the absence 
of strong phases this would result in a cancellation and 
lead to a very small r d . Note that, as in the b — > sj 
case, we cannot predict the sign of cos 5 d . Moreover, 
the values of r d and 5 d are sensitive to the difference 
between the strong phases of the c and u loops, for which 
we do not consider the perturbative result reliable. In 
particular, the matrix element of 0' 2 u in b — > dj may 
have sizable long distance contributions. In any event, 
the short distance calculation predicts that the strong 
phase vanishes for the u loop contribution, while it is 
sizable for the c loop. Therefore, the cancellation in r d is 
unlikely to be effective. The important point is that we 
expect rd ~ ~ r s . 

The crucial difference between the time dependent CP 
asymmetries in B —* f s j and B — > jdl is that in the lat- 
ter case, naively, there are two very strong suppression 
factors. First, considering O-j only, 5 l / d7 is suppressed 
by m d /m b . Second, the phase of the dominant decay 
amplitude, V t bV t * d , cancels the phase of the B° — B° 
mixing amplitude, yielding another strong suppression 
of Sf d7 . Both of these suppressions are fictitious, since 
the V u bV* d (0' 2 — 0' 2u ) contributions lift both the md/mb 
suppression, just like in the b — > case discussed earlier, 
and also the suppression coming from the cancellation of 
the mixing and the decay phases. The leading contribu- 
tion to Sf dl is proportional to r d , which gives a contri- 
bution to Ar with weak phase 7. Using the fact that the 
phase of the mixing amplitude is 2/3, and that of Al is 
/3, we obtain at leading order 



s fdi = -2 



V ubV* d 



" Id 



^cos^sin(/3 + 7 ). (15) 



This result is independent at first order in of the small 
direct CP violation in b — > in the SM. 

A model dependent way to connect the inclusive cal- 
culations with exclusive B — ► K*"f or pry is the Ali— 
Greub model [16] obtained by smearing the inclusive 
rate with a model shape function and integrating over 
E 1 > (m 2 B — lGeV 2 )/(2ms) (i.e., attributing the rate 
to mx < 1 GeV to the K* or p). Restricting the shape 
function parameters to reproduce B(B — > K*j), we ob- 
tain Tk* ~ 0.025, with little sensitivity to the model 
parameters. This is comparable in magnitude but in- 
dependent of the m s jm\, contribution. This model also 
yields r p ~ tk*, which is much larger than md/mb- To 
compare with the "more inclusive" measurement [6] , we 
computed (r s ) mXs <i.8GeV ~ 0.055 using simple shape 
function models. Given the model uncertainties, this fig- 
ure should be taken as a very rough estimate. 



Oi_6 contribute only to Al to all orders in a s at lead- 
ing order in Aqcd /mi (this is viewed as of the same or- 
der as mx* / Ek*)- We identify several types of sublead- 
ing SCET operators that give amplitudes to right-handed 
photons, and yield tk* and r p of order Aqcd/™;>- 

Using SCET, a factorization theorem for heavy-to- 
light form factors has been proven at leading order in 
AQCD/wib [18-20]. There are two contributions to the 
form factors, a nonfactorizablc (or soft or form- factor) 
part and factorizable (or hard scattering) part, which are 
of the same order in Aqqd/to6- 

The exclusive radiative decays considered here were 
analyzed in an expansion in h-QCT>/ m b hi [9, 21-23], and 
recently in SCET [24]. The nonfactorizable part receives 
its dominant contribution from O7 and gives rise to Al 
only. The operator O2 enters at leading order only via 
factorizable contributions. 

When operators O^i are included, the hierarchy of 
the relevant scales gets rather complicated, since m\ > 
m\ ~ AQCDTOb > Aq CD . Ref. [24] assumed that the cc 
loop is dominated by hard loop momenta and can be in- 
tegrated out near the scale mb- Here we adopt a simpler 
approach, by neglecting the charm mass and assuming 
that its effects can be included as a perturbation using 
the formalism of Ref. [25] without encountering singular- 
ities. This assumption is borne out by the explicit onc- 
and two- loop calculations. 

The effective Hamiltonian in Eq. (5) is matched in 
SCETi onto 



G F V t bV t * s e 
HcS ~ V2V ^ 



C(0J) Sn^J^mbPh b v 



(16) 



+ b 1L (cj t )0^ (cji) + b 1R {uJi) (lii) (ui) + C(A 2 ) , 



where the oj's are the usual collinear label momenta. The 
relevant modes are soft quarks and gluons with momenta 
k s ~ A and two types of collinear quarks and gluons 
along n and n (including charm, which can be soft or 
collinear) . Note that the approach of Ref. [24] , treating 
m c ~ m,b, will likely require additional modes. We take 
the photon momentum as q p = E^n^, the collinear s 
quark to move along n^, and A p denotes the transverse 
photon field. The operator in the first line in Eq. (16) 
occurs at leading order in the expansion parameter, A = 
y/A/mb, and its Wilson coefficient is dominated by CV, 
c(ui) = Cr + 0[a s (mb)]- The operators in the second line, 



n ■ V 



P R b Vl 



-^igfa j^Pubv, (17) 



III. EXCLUSIVE B — > K*f AND B -> pj 

We consider next the photon polarization in exclusive 
B — > K*j and pj decays using SCET [17]. We prove that 



are the only SCETi operators suppressed by A that cou- 
ple to a transverse photon and are allowed by power 
counting and s chirality. Here igB% = [n ■ iD c , i-DJfJ is 
the collinear gluon field tensor; for the remaining nota- 
tions see Ref. [26]. The operators 0( 1L ' and 0( lfl ) couple 
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only to jl,r, respectively. Their Wilson coefficients are such as the C(A 2 ) operator 



&ii,(wi,W2) = C 7 + C 2 /3 + 0[C3-6,a s (m b )}, 
6i fl (wi,W2) = -C , 2/3 + 0[C , 3-6,a.(m 6 )]. (18) 

Although the operators in the first and second lines 
of Eq. (16) are of different orders in A, after match- 
ing onto SCETn, they contribute at the same order in 
Aqcd/w6 to the B — > K*-f amplitude. The leading order 
c(w) term gives the nonfactorizable contribution, which 
only contributes to Al, while C>( 1L ) and 0^ 1R ^ give the 
factorizable contributions through time-ordered products 
with the ultrasoft-collinear Lagrangian, 

4, [27]. After 
matching onto SCETn one finds, schematically 

J^xT{0^ L \wi),iC^(x)} - Jdk + J ± (4\k + ) 
® [(gK n ) fc+ fM ± 7±i , ii(i r n + 6»)] (V^ft7±«n,a,i),(19) 
r d 4 a;T{0( 1,l )(w < ),i4i ) (a:)} - J dk+ J\\(u^ , k+) 
® [(gy„) fe+ ^ ± P R (y n t 6„)] (s^Pl^), (20) 

where J7jii are jet functions that have expansions in 
a s (^/AQCDW(,). The operator in Eq. (19) contributes 
to B — > i^*7 at leading order, 



V2tt 2 2 



c(m B )C. 



+ 



m B 



dxdtdk i 



x &il(tob(1 - t),m B t) J±(x,t,k + )(j)K,{x) (j)^{k + ) ,(21) 

where all nonperturbative matrix elements are defined 
as in [26]. The operator in Eq. (20) only gives rise to 
longitudinal K* and therefore does not contribute to 
B — > if*7. This proves that (9 2 contributes at lead- 
ing order in AQCT>/m b to B — > if*7L via Eq. (19), but 
its contribution to £> — > K*jr vanishes at this order. 
(Interestingly, it would contribute to B* — > K^*^r at 
leading order.) The same proof also holds for the other 
four-quark operators, 0±-6- 

This result agrees with the 0(a s ) computation [21, 22], 
and extends it to all orders in a s . This is important, 
since cc s (/x) at the hard-collinear scale /i 2 ~ Am b may or 
may not be perturbative [18, 28, 29]. The suppression of 
Ar can be understood from a simple hclicity argument. 
Since the s quark is left-handed, the K* cannot be right- 
handed, unless additional right-handed gluons end up in 
the final meson. However, the contributions of higher 
Fock states are power suppressed. 

The B — > K*jr amplitude does arise at subleading or- 
der in Aqcd/^6- There are several sources of such cor- 
rections. For example, (i) time ordered products of 0^ 1R ^ 
with the subleading collinear Lagrangian which 
lead to factorizable contributions similar to Eq. (19) con- 
taining an explicit factor a s (^/AqcD m b); (ii) higher or- 
der terms in the SCETi effective Lagrangian in Eq. (16), 



(2R \uj) 



dx+ 
~2^~ 



dk _ e -ix+k-/2 i 



2£LA;_ 



x ii n „[Yiip™Y n ](^) j ± P R b v , (22) 

where the gauge invariant operator \Y^ip ns Y n }{x) con- 
tains Wilson lines of the ultrasoft fields in the n direc- 
tion. 0( 2fl ) is obtained by matching the graph in Fig. 1 
with one ultrasoft gluon, and its Wilson coefficient is re- 
lated to bin in Eq- (18) by reparameterization invari- 
ance [30]. The C 2 /3 terms in Eq. (18) correspond to 
m c = 0, when k(z) = 1/2. In this case the fc_ and x + 
integrals in Eq. (22) can be performed trivially. How- 
ever, in Eq. (22) we kept m c ^ to exhibit the analogy 
with the nonperturbative corrections to inclusive b — ► sj 
decay [31, 32]. For m c ^ 0, additional nonlocality is in- 
troduced by k(z) = 1/2 — 2arctan[y / z/(4— z)] 2 /z with 
z = 2E J k-/ml [12]. [For m c ^ 0, the C 2 /3 terms in 
Eq. (18) should be multiplied by 2K{2E 1 Lo 2 /m 2 c )] Ex- 
panding 0( 2K ) in x + reproduces the series of operators 
proportional to {m b AQcv I m V) n in Eq. (5) in [32]. Time- 
ordered products of C>( 2fl ) with the ultrasoft-collinear La- 

(2) 

grangian, Q , give 0(Aqcd/?™!>) soft contributions to 
Ar, with no a s suppression compared to Al. 

A complete study of these subleading contributions is 
rather involved and we leave it for future work. It seems 
unlikely to us that a cancellation could result in a sup- 
pression of r K » and r p to order Aq CD /to 2 . This leads to 
the dimensional estimate 



r K , 



1 C 2 A 



QCD 



3 CV m b 



0.1. 



(23) 



This effect dominates over the m s piece of O7 for tk* 
The estimate for r p is more involved because of the con- 
tributions with different weak phases, 



r K * 



1 



Vu h v; d 

Vt b v; d 



loop 



hit 



CwA 47T • 



A 



QCD 

m b 



(24) 

The Cioop term comes from the non-cancellation of the c 
and u loops, and we expect it to have a numerically large 
coefficient. The Cwa term arises from weak annihilation, 
whose contribution to Ar at order Aqcd / m b vanishes [9] . 
The latter contributes significantly only in B ± , while in 
B° it is color suppressed. Thus, we expect that the SM 
prediction for S pl is not much smaller than it is for Sk*j- 
For higher mass one-body hadronic final states, Ar 
still vanishes at leading order, but the suppression by 
mx /Ex is expected to be less effective as mx increases 
(although there is no evidence for this in the B —> DX 
data [33]). Thus, the SM value of r is expected to de- 
pend on the final state. For high-mass and multi-body 
final states Ar may arise, formally, at leading order in 
Aqcd/™&- For example, B — ► B*tt^ so{ ^ followed by 
B* — > K^jr can give rise to B° — ► K s tt°jr with 
m K s ir a ~ \/Aqcd?ti&, without a Aqcd/tJ() suppression. 
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Therefore, averaging the results of B — > K*j [3, 4] with 
B — > Kstt°j [6] is not free from theoretical uncertainties. 

IV. DISCUSSION AND SUMMARY 

We studied the standard model prediction for the _B — > 
X s ,djR decay amplitude with a right-handed photon, An, 
compared to the amplitude with a left-handed photon. 
Considering only O7, their ratio is r Si d = m s ^/m}>, which 
reproduces the often-quoted prediction independent of 
the hadronic final state. However, including the other 
operators, Oi^7, An becomes much larger than this naive 
estimate, and hadronic physics gives rise to sizable un- 
certainties. The time dependent CP asymmetries also 
become sensitive to the strong phase. 

In inclusive B — ► X s -f and X^n decays, An is only 
suppressed by g s /(A-n). We calculated r inclusively, and 
found it to be of order 0.1 depending on the cut on the 
photon energy [see Eq. (8)]. While this calculation is reli- 
able, it cannot be easily compared to data. If one restricts 
the hadronic final state, such as in the measurement of 
the time dependent CP asymmetry Sk s h°~( with a range 
of Kstt invariant masses, then that can no longer be 
calculated using inclusive methods. Still, our results in- 
dicate that it would be hard to argue that a measurement 
of S Ks1T o 7 ~ 0.1 cannot be due to SM physics. 

In exclusive B — > K*j and p-f decays, we proved using 
SCET that An vanishes at lowest order in the Aqcd / rrib 
expansion to all orders in a s . The leading contribution to 
An is formally of order Aqcd/w;,, but numerically it is 
enhanced by C 2 /C 7 [see Eqs. (23) and (24)]. The result 
depends on unknown hadronic matrix elements, which 



give rise to a sizable uncertainty in the SM prediction. 

Both the inclusive and the exclusive calculations pre- 
dict that r s and are comparable to each other in 
the SM. Thus, we also expect the time dependent CP 
asymmetries in B° — > /<j7 to be comparable to those in 
B° — > / s 7 (except for a modest suppression by IV^/V^I), 
in contrast to what has been widely believed. We esti- 
mate these asymmetries in decays such as B — > AT*7, 
P7, K s n°j, and ir + n~j to be of order 0.1, with large 
uncertainties. 

To conclude, our main result is that the standard 
model prediction for Ar/Al is of order 0.1, with sizable 
uncertainty. A measurement of it much above this level 
would indicate new physics. To draw conclusions from a 
smaller value would require a more complete analysis and 
knowledge of hadronic matrix elements including strong 
phases. More effort in this direction would be welcome. 
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